Panmixia in the American eel extends to its tropical range of distribution: Biological implications and policymaking challenges.

The American eel (Anguilla rostrata) has long been regarded as a panmictic fish and has been confirmed as such in the northern part of its range. In this paper, we tested for the first time whether panmixia extends to the tropical range of the species. To do so, we first assembled a reference genome (975 Mbp, 19 chromosomes) combining long (PacBio and Nanopore and short (Illumina paired-end) reads technologies to support both this study and future research. To test for population structure, we estimated genotype likelihoods from low-coverage whole-genome sequencing of 460 American eels, collected at 21 sampling sites (in seven geographic regions) ranging from Canada to Trinidad and Tobago. We estimated genetic distance between regions, performed ADMIXTURE-like clustering analysis and multivariate analysis, and found no evidence of population structure, thus confirming that panmixia extends to the tropical range of the species. In addition, two genomic regions with putative inversions were observed, both geographically widespread and present at similar frequencies in all regions. We discuss the implications of lack of genetic population structure for the species. Our results are key for the future genomic research in the American eel and the implementation of conservation measures throughout its geographic range. Additionally, our results can be applied to fisheries management and aquaculture of the species.

A qPCR-based method to detect the eel parasitic nematode Anguillicola crassus in intermediate and final hosts.

Being able to systematically detect parasitic infection, even when no visual signs of infection are present, is crucial to the establishment of accurate conservation policies. The nematode Anguillicola crassus infects the swimbladder of anguillid species and is a potential threat for eel populations. In North America, naïve hosts such as the American eel Anguilla rostrata are affected by this infection. The accidental introduction of A. crassus following restocking programs may contribute to the actual decline of the American eel in Canada. We present a quantitative real time PCR-based method to detect A. crassus infection in final and intermediate hosts. We tested two protocols on samples from different geographical origins in Canada: 1) a general detection of A. crassus DNA in pools of young final hosts (glass eels) or crustacean intermediate hosts 2) a detection at the individual scale by analyzing swim bladders from elvers, or from adult yellow and silver eels. The DNA of A. crassus was detected in one pool of zooplankton (intermediate host) from the Richelieu River (Montérégie-Québec), as well as in individual swim bladders of 13 elvers from Grande and Petite Trinité rivers (Côte-Nord-Québec). We suggest that our qPCR approach could be used in a quantitative way to estimate the parasitic burden in individual swim bladders of elvers. Our method, which goes beyond most of previous developed protocols that restricted the diagnosis of A. crassus to the moment when it was fully established in its final host, should help to detect early A. crassus infection in nature.

Long-distance migration is a major factor driving local adaptation at continental scale in Coho salmon.

Inferring the genomic basis of local adaptation is a long-standing goal of evolutionary biology. Beyond its fundamental evolutionary implications, such knowledge can guide conservation decisions for populations of conservation and management concern. Here, we investigated the genomic basis of local adaptation in the Coho salmon (Oncorhynchus kisutch) across its entire North American range. We hypothesized that extensive spatial variation in environmental conditions and the species' homing behaviour may promote the establishment of local adaptation. We genotyped 7829 individuals representing 217 sampling locations at more than 100,000 high-quality RADseq loci to investigate how recombination might affect the detection of loci putatively under selection and took advantage of the precise description of the demographic history of the species from our previous work to draw accurate population genomic inferences about local adaptation. The results indicated that genetic differentiation scans and genetic-environment association analyses were both significantly affected by variation in recombination rate as low recombination regions displayed an increased number of outliers. By taking these confounding factors into consideration, we revealed that migration distance was the primary selective factor driving local adaptation and partial parallel divergence among distant populations. Moreover, we identified several candidate single nucleotide polymorphisms associated with long-distance migration and altitude including a gene known to be involved in adaptation to altitude in other species. The evolutionary implications of our findings are discussed along with conservation applications.

Re-evaluating Coho salmon (Oncorhynchus kisutch) conservation units in Canada using genomic data.

Conservation units (CUs) are important tools for supporting the implementation of standardized management practices for exploited species. Following the adoption of the Wild Salmon Policy in Canada, CUs were defined for Pacific salmon based on characteristics related to ecotype, life history and genetic variation using microsatellite markers as indirect measures of local adaptation. Genomic data sets have the potential to improve the definition of CUs by reducing variance around estimates of population genetic parameters, thereby increasing the power to detect more subtle patterns of population genetic structure and by providing an opportunity to incorporate adaptive information more directly with the identification of variants putatively under selection. We used one of the largest genomic data sets recently published for a nonmodel species, comprising 5662 individual Coho salmon (Oncorhynchus kisutch) from 149 sampling locations and a total of 24,542 high-quality SNPs obtained using genotyping-by-sequencing and mapped to the Coho salmon reference genome to (1) evaluate the current delineation of CUs for Coho in Canada and (2) compare patterns of population structure observed using neutral and outlier loci from genotype-environment association analyses to determine whether separate CUs that capture adaptive diversity are needed. Our results reflected CU boundaries on the whole, with the majority of sampling locations managed in the same CU clustering together within genetic groups. However, additional groups that are not currently represented by CUs were also uncovered. We observed considerable overlap in the genetic clusters identified using neutral or candidate loci, indicating a general congruence in patterns of genetic variation driven by local adaptation and gene flow in this species. Consequently, we suggest that the current CU boundaries for Coho salmon are largely well-suited for meeting the Canadian Wild Salmon Policy's objective of defining biologically distinct groups, but we highlight specific areas where CU boundaries may be refined.

Adaptive and maladaptive genetic diversity in small populations: Insights from the Brook Charr (Salvelinus fontinalis) case study.

Investigating the relative importance of neutral versus selective processes governing the accumulation of genetic variants is a key goal in both evolutionary and conservation biology. This is particularly true in the context of small populations, where genetic drift can counteract the effect of selection. Using Brook Charr (Salvelinus fontinalis) from Québec, Canada, as a case study, we investigated the importance of demographic versus selective processes governing the accumulation of both adaptive and maladaptive mutations in closed versus open and connected populations to assess gene flow effect. This was achieved by using 14,779 high-quality filtered SNPs genotyped among 1,416 fish representing 50 populations from three life history types: lacustrine (closed populations), riverine and anadromous (connected populations). Using the PROVEAN algorithm, we observed a considerable accumulation of putative deleterious mutations across populations. The absence of correlation between the occurrence of putatively beneficial or deleterious mutations and local recombination rate supports the hypothesis that genetic drift might be the main driver of the accumulation of such variants. However, despite a lower genetic diversity observed in lacustrine than in riverine or anadromous populations, lacustrine populations do not exhibit more deleterious mutations than the two other history types, suggesting that the negative effect of genetic drift in lacustrine populations may be mitigated by that of relaxed purifying selection. Moreover, we also identified genomic regions associated with anadromy, as well as an overrepresentation of transposable elements associated with variation in environmental variables, thus supporting the importance of transposable elements in adaptation.

Combined effects of temperature changes and metal contamination at different levels of biological organization in yellow perch.

In this study, we measured the effects of temperature (9°C, 20°C, and 28°C), metal contamination (cadmium and nickel) and their interaction on yellow perch (Perca flavescens) using liver enzymatic and transcriptomic endpoints and biometric indices. Kidney metal concentrations increased with a rise of temperature. The biometric indices analysed (Fulton condition factor, pyloric cæca, hepatosomatic and gonadosomatic indices) generally decreased with an increase of temperature but not with metal contamination. At the enzymatic level, the activity of superoxide dismutase (SOD), involved in antioxidant response, was affected by both temperature and metal contamination, whereas the activity of glucose-6-phosphate dehydrogenase (G6PDH), involved in energy accumulation but also in antioxidant response, was only affected by metal exposure. The response of perch to the stressors at the transcriptional level differed from the metabolic response. In particular, the transcription level of the cco and g6pdh genes sharply decreased with increasing temperature, while the activities of the corresponding enzymes remained stable. The normal response of the transcription level of the apoptotic gene (diablo) to heat stress was also altered in metal-contaminated fish. The combination of metal and temperature stresses also modified the response of antioxidant metabolism induced by these stressors individually. This study contributes to a better understanding of the influences of natural stressors like temperature on biomarkers commonly used in ecotoxicological studies and will facilitate their interpretation in the context of multiple stressors characteristic of field situations.

Transcriptional response of yellow perch to changes in ambient metal concentrations-A reciprocal field transplantation experiment.

Recent local adaptation to pollution has been evidenced in several organisms inhabiting environments heavily contaminated by metals. Nevertheless, the molecular mechanisms underlying adaptation to high metal concentrations are poorly understood, especially in fishes. Yellow perch (Perca flavescens) populations from lakes in the mining area of Rouyn-Noranda (QC, Canada) have been faced with metal contamination for about 90 years. Here, we examine gene transcription patterns of fish reciprocally transplanted between a reference and a metal-contaminated lake and also fish caged in their native lake. After four weeks, 111 genes were differentially transcribed in metal-naïve fish transferred to the metal-contaminated lake, revealing a plastic response to metal exposure. Genes involved in the citric cycle and beta-oxidation pathways were under-transcribed, suggesting a potential strategy to mitigate the effects of metal stress by reducing energy turnover. However, metal-contaminated fish transplanted to the reference lake did not show any transcriptomic response, indicating a reduced plastic response capability to sudden reduction in metal concentrations. Moreover, the transcription of other genes, especially ones involved in energy metabolism, was affected by caging. Overall, our results highlight environmental stress response mechanisms in yellow perch at the transcriptomic level and support a rapid adaptive response to metal exposure through genetic assimilation.

Dietary sodium protects fish against copper-induced olfactory impairment.

Exposure to low concentrations of copper impairs olfaction in fish. To determine the transcriptional changes in the olfactory epithelium induced by copper exposure, wild yellow perch (Perca flavescens) were exposed to 20 µg/L of copper for 3 and 24h. A novel yellow perch microarray with 1000 candidate genes was used to measure differential gene transcription in the olfactory epithelium. While three hours of exposure to copper changed the transcription of only one gene, the transcriptions of 70 genes were changed after 24h of exposure to copper. Real-time PCR was utilized to determine the effect of exposure duration on two specific genes of interest, two sub-units of Na/K-ATPase. At 24 and 48 h, Na/K-ATPase transcription was down-regulated by copper at olfactory rosettes. As copper-induced impairment of Na/K-ATPase activity in gills can be ameliorated by increased dietary sodium, rainbow trout (Oncorhynchus mykiss) were used to determine if elevated dietary sodium was also protective against copper-induced olfactory impairment. Measurement of the olfactory response of rainbow trout using electro-olfactography demonstrated that sodium was protective of copper-induced olfactory dysfunction. This work demonstrates that the transcriptions of both subunits of Na/K-ATPase in the olfactory epithelium of fish are affected by Cu exposure, and that dietary Na protects against Cu-induced olfactory dysfunction.

A multi-level biological approach to evaluate impacts of a major municipal effluent in wild St. Lawrence River yellow perch (Perca flavescens).

The development of integrated ecotoxicological approaches is of great interest in the investigation of global concerns such as impacts of municipal wastewater effluents on aquatic ecosystems. The objective of this study was to investigate the effects of a major wastewater municipal effluent on fish using a multi-level biological approach, from gene transcription and enzyme activities to histological changes. Yellow perch (Perca flavescens) were selected based on their wide distribution, their commercial and recreational importance, and the availability of a customized microarray. Yellow perch were sampled upstream of a major municipal wastewater treatment plant (WWTP) and 4 km and 10 km downstream from its point of discharge in the St. Lawrence River (Quebec, Canada). Concentrations of perfluoroalkyl substances (PFASs), polybrominated diphenyl ethers (PBDEs) and metals/trace elements in whole body homogenates were comparable to those from other industrialized regions of the world. Genomic results indicated that the transcription level of 177 genes was significantly different (p<0.024) between exposed and non-exposed fish. Among these genes, 38 were found to be differentially transcribed at both downstream sites. Impacted genes were associated with biological processes and molecular functions such as immunity, detoxification, lipid metabolism/energy homeostasis (e.g., peroxisome proliferation), and retinol metabolism suggesting impact of WWTP on these systems. Moreover, antioxidant enzyme activities were more elevated in perch collected at the 4 km site. Biomarkers of lipid metabolism, biosynthetic activity, and aerobic capacities were significantly lower (p<0.05) in fish residing near the outfall of the effluent. Histological examination of the liver indicated no differences between sites. Correlations between PFAS, PBDE, and metal/trace element tissue concentrations and markers of peroxisomal proliferation, oxidative stress, and retinoid metabolism were found at the gene and cellular levels. Present results suggest that relating transcriptomic analyses to phenotypic responses is important to better understand impacts of environmental contamination on wild fish populations.

Chemosensory mediated behaviors and gene transcription profiles in wild yellow perch (Perca flavescens) from metal contaminated lakes.

The olfactory system of fish is sensitive to the toxic effects of low concentrations of contaminants. To investigate the effects of long-term metal exposure on olfaction in wild yellow perch (Perca flavescens), fish from one clean (Geneva Lake) and two metal-contaminated lakes (Ramsey and Hannah lakes) were collected in and around the metal-mining district of Sudbury, ON. Two different techniques were used to measure the effects of exposure to environmental contamination: (i) behavioral responses were recorded in response to conspecific skin extract and (ii) gene transcription differences in olfactory rosettes were characterized using a novel, 1000-candidate gene yellow perch microarray. Behavioral assays performed on fish from the clean lake demonstrated avoidance of a conspecific skin extract, while fish from metal contaminated lakes showed no avoidance response. A total of 109 out of the 1000 genes were differentially transcribed among the lakes. Most of the differentially transcribed genes were between the two metal contaminated lakes relative to either of the contaminated lakes and the reference lake. No genes were differentially expressed between Geneva Lake (clean) and Hannah Lake (metal contaminated). These results demonstrated that even though the different populations of fish from both Hannah and Ramey lakes were affected at the behavioral level, the impairment of olfaction was not measurable using gene transcriptional changes in olfactory rosettes.

How does exposure to nickel and cadmium affect the transcriptome of yellow perch (Perca flavescens)--results from a 1000 candidate-gene microarray.

The molecular mechanisms underlying nickel (Ni) and cadmium (Cd) toxicity and their specific effects on fish are poorly understood. Documenting gene transcription profiles offers a powerful approach toward identifying the molecular mechanisms affected by these metals and to discover biomarkers of their toxicity. However, confounding environmental factors can complicate the interpretation of the results and the detection of biomarkers for fish captured in their natural environment. In the present study, a 1000 candidate-gene microarray, developed from a previous RNA-seq study on a subset of individual fish from contrasting level of metal contamination, was used to investigate the transcriptional response to metal (Ni and Cd) and non metal (temperature, oxygen, and diet) stressors in yellow perch (Perca flavescens). Specifically, we aimed at (1) identifying transcriptional signatures specific to Ni and Cd exposure, (2) investigating the mechanisms of their toxicity, and (3) developing a predictive tool to identify the sublethal effects of Ni and Cd contaminants in fish sampled from natural environments. A total of 475 genes displayed significantly different transcription levels when temperature varied while 287 and 176 genes were differentially transcribed at different concentrations of Ni and Cd, respectively. These metals were found to mainly affect the transcription level of genes involved in iron metabolism, transcriptional and translational processes, vitamin metabolism, blood coagulation, and calcium transport. In addition, a linear discriminant analysis (LDA) made using gene transcription levels yielded 94% correctly reassigned samples regarding their level of metal contamination, which indicates the potential of the microarray to detect perch response to Cd or Ni effects.

The transcriptional landscape of cross-specific hybrids and its possible link with growth in brook charr (Salvelinus fontinalis Mitchill).

The genetic mechanisms underlying hybridization are poorly understood despite their potentially important roles in speciation processes, adaptative evolution, and agronomical innovation. In this study, transcription profiles were compared among three populations of brook charr and their hybrids using microarrays to assess the influence of hybrid origin on modes of transcription regulation inheritance and on the mechanisms underlying growth. We found that twice as many transcripts were differently expressed between the domestic population and the two wild populations (Rupert and Laval) than between wild ones, despite their deeper genetic distance. This could reflect the consequence of artificial selection during domestication. We detected that hybrids exhibited strikingly different patterns of mode of transcription regulation, being mostly additive (94%) for domestic × Rupert, and nonadditive for Laval × domestic (45.7%) and Rupert × Laval hybrids (37.5%). Both heterosis and outbreeding depression for growth were observed among the crosses. Our results indicated that prevalence of dominance in transcription regulation seems related to growth heterosis, while prevalence of transgressive transcription regulation may be more related to outbreeding depression. Our study clearly shows, for the first time in vertebrates, that the consequences of hybridization on both the transcriptome level and the phenotype are highly dependent on the specific genetic architectures of crossed populations and therefore hardly predictable.